MPP-2009-33 



O 
O 



Multiple Spectral Splits of Supernova Neutrinos 

Basudcb Dasgupta/ Amol Dighc,^ Gcorg G. RafFclt/ and Alcxci Yu. Smirnov^'^ 

' Max-Planck-Institut fiir Physik (Werner-Heisenberg-Institut), Fdhringer Ring 6, 80805 Munchen, Germany 
^Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400005, India 
'^Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34014 Trieste, Italy 
'^Institute for Nuclear Research, Russian Academy of Sciences, 117312 Moscow, Russia 
(Dated: 19 April 2009, revised 17 June 2009) 

Collective oscillations of supernova neutrinos swap the spectra f^^iE) and fv^{E) with those of 
another flavor in certain energy intervals bounded by sharp spectral splits. This phenomenon is far 
more general than previously appreciated: typically one finds one or more swaps and accompanying 
splits in the and v channels for both inverted and normal neutrino mass hierarchies. Depending 
on an instability condition, swaps develop around spectral crossings (energies where /^^ = fv^ , 
fv^ = fa^ as well as i5 oo where all fluxes vanish), and the widths of swaps are determined by 
the spectra and fluxes. Wash-out by multi-angle decoherence varies across the spectrum and splits 
can survive as sharp spectral features. 

PACS numbers: 14.60.Pq, 97.60.Bw 
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Introduction. — The neutrino flux from a core-collapse 
supernova (SN) is a powerful probe of particle physics 
and astrophysics [l[ . SN neutrinos interact not only with 
the stellar medium, producing the Mikheyev-Smirnov- 
Wolfenstein (MSW) flavor conversion, but also with other 
neutrinos and antineutrinos. The latter interactions 
modify the flavor evolution in a non-linear fashion and 
give rise to collective forms of oscillation B, S 3, B, 3 , 
a subject of intense recent investi gati on 
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The most important observational consequence of the 
collective effects is an exchange of the Ve (i^e) spec- 
trum with the Vx {vx) spectrum in certain energy inter- 
vals. We call such a flavor exchange a "swap" , whereas 
"splits" are sharp boundary features at the edges of each 
swap interval. Spectral splits may become observable in 
the high-statistics neutrino signal from the next galac- 
tic SN, leading to valuable clues about the underlying 

physics [H,!!!,!!!]. 



The well-understood "classic swap" covers the entire 
v spectrum and that of v above an energy fixed by the 
approximate conservation of the Ve deleptonization flux 
la . 17, In this paper we show that spectral swaps 



and concomitant splits are more ubiquitous than has 
been appreciated in the past. One example is the puz- 
zling low-energy split in the 9 spectrum that was noted 



for the inverted neutrino mass hierarchy [20|, |2l| . How- 
ever, with flavor spectra typical for SN neutrinos one 
should expect multiple splits in either hierarchy. 

We focus on neutrino-neutrino interactions alone and 
study two-flavor oscillations driven by the atmospheric 
mass difference and 1-3 mixing. As has been established 
before [lj| , the usual matter effect in the region of collec- 
tive oscillations (up to a few 100 km) can be accounted for 
by choosing a small (matter suppressed) effective mixing 
angle which we take to be Ocs = 10~^. MSW conver- 
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FIG. 1: SN neutrino spectra before (dashed lines) and after 
(solid lines) collective oscillations, but before possible MSW 
conversions. The panels are for v and v, each time for IH and 
NH. Red lines e-flavor, blue i-flavor. Shaded regions mark 
swap intervals. 



sions occur typically at larger distances. Their effects 



then factorize and can be included separately [24 1. 

Spectral crossings and spectral swaps. — Consider first 
the SN cooling phase where plausible choices are [s^l 
F^^ : % : F^^ =_0.85 : 0.75 : _1.00 for the neutrino 
fluxes, eI^ = 12, E^^ = 15 and E^^ = E^^ =18 MeV 
for the average energies, and fu{E) cx E^e~^^^^ for the 
spectral shape. Based on the single-angle approximation 
for neutrino propagation [ill. [l3. [20I. |28|| Fig. [T] shows the 
flavor spectra before and after collective oscillations. For 
the inverted mass hierarchy (IH) we find a swap for both 
v and V and thus a total of four splits. For the normal 
hierarchy (NH) the swaps extend to infinite E, providing 
one split in the ly and P spectrum each. 

Flavor oscillations leave f^^ {E) -\- f^^ {E) unchanged. 
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FIG. 2: Initial spectrum g^j defined in Eq. ([ij for the cooling- 
phase example of Fig. [1] Swap factor after collective oscil- 
lations for IH and NH as indicated. Single-an gle (solid) and 
multi-angle (dashed). For animations see Ref. [33| . 



so Af,y{E) = fuA^) ~ fv^i^); s-iid similarly for v, con- 
tains all relevant information. It proves crucial for under- 
standing the multisplit phenomenon (i) to use —Afj^{E) 
for 9 (the "flavor isospin" construction [ll[), and (ii) to 
merge the v and 9 spectra to a single continuum in terms 
of the variable lu = ±|Am^/2£^|, the vacuum oscilla- 
tion frequency, that is positive for v and negative for 
£/. Here Am^ = m| — mf with vi w and ;/2 ~ i^x, 
so Am^ > stands for NH and Am^ < for IH. With 
~ \A.vr? wc thus define our spectrum as 
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(1) 



where the common factor comes from the E ^ lo trans- 
formation. In SN we have f^^ (E) = /p^ [E) and a net 

dcleptonization flux, so ^J^Aiog^j > 0. Notice that 
g^} > means excess z/g for to > and excess for uj < 0. 

In Fig.[2]we show g{uj) corresponding to the example of 
Fig. [H The effect of collective oscillations is represented 
by a "swap factor" S^, giving the final spectrum as 



^ final 



q ini 



tial 



(2) 



The "spectral crossings" where g{Lo) ~ are indicated 
by vertical dotted lines. In particular, our construction 
provides a ly-D crossing at w = (infinite E), playing 
the same role as a "true flavor crossing" where /jy^ (E) = 
fi,^ (E) . The v-D crossing exists even if the source emits 
only i/g and Pg. 



FIG. 3: Same as Fig. [2] for the accretion-phase example of 
Refs. [13, [2l[. Thin solid line: swap factor for e^s = 10~^° 
instead of 10~^. There is no resolved swap for NH. 



Figure [2] illustrates a general point: for IH every g^ 
crossing with a positive slope is an "unstable crossing" , 
i.e. it induces a swap over a finite to range with splits 
that become sharp in the adiabatic limit. For NH every 
g^^ crossing with a negative slope is an unstable cross- 
ing. (The equations of motion actually reveal that the 
dynamics for a given spectrum g^ in IH is equivalent to 
a spectrum —gi^ in NH, thus there is no fundamental 
difference between the two hierarchies.) 

The usual assumptions about SN neutrinos lead to a 
triple-crossed spectrum and thus to one or two swaps 
(two or four splits), depending on the hierarchy. An ap- 
parent counter example, representing the SN accretion 
phase, is E^^ ^ 10, Ep^ = 15 and E^^ = = 24 MeV 
and F^^ : fI^ : F^^ =^2.4 : 1.6 : l.o'ji^, [iJ. The cor- 
responding iu spectrum (Fig. [3]) has three crossings, yet 
only one swap is observed for IH and none for NH. Below 
we argue that the narrow spacing of the crossings implies 
that (i) in the adiabatic limit the missing swap is expo- 
nentially narrow, and (ii) it is suppressed by adiabaticity 
violation. Therefore, a narrowly spaced triple crossing 
can superficially act like a single one. 

Single-crossed system. — Many of our general observa- 
tions arc borne out from a generic single-crossed example 
for IH (Fig. [3]). The equations of motion (EOMs) are 



{loB + ^P) X 



(3) 



where B is a unit vector (0,0, 1) in flavor space, P^ is 
the polarization-vector spectrum that covers v for uj > 
and D for ui < 0, and P = dwP^j. The coupling /j, 
depends on Gp, the local v and D densities, and the nor- 
malization of g^. Initially P^^ = (0,0, g^^), for a vanish- 
ingly small mixing angle. The crossing is unstable for IH 
because a swap would decrease the energy / dw wB • P;^ . 
(For NH, the wB term in Eq. ^ appears with a nega- 
tive sign, so J dujLo'B ■ is already minimal. However, 
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FIG. 4: Basic case of a spectrum producing a swap. The ini- 
tial spectrum (dashed) has an unstable crossing at Uc, whereas 
the final spectrum (solid) has splits at uJi and uj2- 



a spectrum —g^^ would have an unstable crossing.) We 
seek the final spectrum for this schematic case, when /i 
slowly decreases from some initial value to 0. 

The EOMs imply that B • P = j dui g^j is conserved, so 
we can "rotate away" a common precession of the entire 
system with frequency ^B-P in the same way as removing 
the usual matter term 1271 . The EOMs become 



(t^B + /iP_L) X P, 



(4) 



where Pj^ is the component of P transverse to B. 

For B • P 7^ there is a critical value /isync such that for 
Mo > /Usync the system initially performs synchronized os- 
cillations with ^Pj^ =0. As /i decreases, each P^; follows 
its Hamiltonian wB -I- //P±, performing a "pure preces- 
sion," and all P^ with lo < tJspiit end anti-aligncd with 
B, the others aligned, leading to the classic split fixed by 
B-P conservation [16, 17, 18J. For B-P = and Up oo, 
the system evolves like the flavor pendulum , ending 
with the entire spectrum swapped. 

A qualitatively different behavior occurs for fiQ < 
Msync, or any finite if B • P = 0. Now fJ.P± performs 
a fast non-adiabatic motion with large amplitude, caus- 
ing a collective oscillation of the entire spectrum. This 
generic behavior is borne out analytically for a spectrum 
that is antisymmetric in u relative to Wc, leading to a 
"pure nutation" in terms of a variable (Mt). Without 
loss of generality we use Wc = and find [30, [HI 




— — . ujK^2{c - c„0 I (5) 

where s = simp, c = cos(y5, Cm = cos(y9,„ax, and ip{t) 
obeys ^(p'^ = K^^cosip — cos (/Smax)- Differentiation pro- 
vides ip ~ — sin(/3, the equation for an anharmonic os- 
cillator. Equation ([5]) solves the EOM (0]) if the nutation 
amplitude </3max and frequency k fulfill the consistency 
relation 



dw 



- 2uj'^K^C^ -t- 



1 



(6) 



For a "symmetric box" where giio) = l/Aw for < w < 
iAw and g^^ = -g^ we find kq = l^uj j \J I ^^^^ - 1. 

The initial state when /i = /ip is described by (yS = 
<i5max ~ TT, implying c,„ = — 1 and so the amplitude of 
Pf, in Eq. (O is Kq/(w^ -f Kq). The modulation of g^ 
thus has a resonance with a maximum at w = 0. As /i 
decreases adiabatically from /zq to 0, the r.h.s. of Eq. ^ 
becomes singular and so the l.h.s. must become singular 
as well. This is only possible for Cm ^ +1, i.e., the fi- 
nal maximum nutation amplitude is i^Smax = 0. At the 
same time, the frequency k decreases from to some 
non-zero value Kg- With c,n = +1, Eq. ([5]) provides 
~ [llP' — Kg)/|w^ — Kg I = ±1, and thus there is a swap 
in the range of frequencies < . In other words, the 
final K plays the role of the split frequency Wsput ~ ±Ks, 
which is related to kq through an adiabatic invariant. 

Because kq plays the dual role of the approximate 
swap width and an inverse evolutionary time scale, an 
exponentially small kq implies that the corresponding 
swap disappears (i) by not being numerically resolved 
and (ii) by adiabaticity violation. 

A non-zero mixing (0 < ^eff <C 1) is required to trig- 
ger the initial motion, but also causes a delay until the 
system relaxes to the adiabatic solution. Therefore, the 
effective relevant for swap formation is reduced and 
the swap width is smaller for a smaller 6*69. This sub- 
tle effect can become visible when a split falls into a tail 
of guj where the flux is small. A case in point is the 
left split in Fig. [3] that shifts to the right for a smaller 
0oft (thin- line) , explaining the dependence on the matter 
profile (that modifies 6'cft) noted in Refs. 2^, [in. Similar 
remarks pertain to the left-most IH split in Fig. [2l 

Multiple crossings. — For a spectrum with several cross- 
ings the evolution is more complicated and depends on 
various features, notably the separation of the crossings 
and ijlq. For instance, a multi-crossed system can be con- 
structed from two clones of Fig. S] with a large u range 
of empty modes between them. If is small, each indi- 
vidual block shows the behavior discussed above. Even 
though the entire spectrum feels the same /iP^(t), only 
the Fourier components close to a given uj strongly affect 
, so different regions of the spectrum can evolve almost 
independently of each other. If /io > ^sync of a single 
clone, we get the "co-rotating plane solution" and clas- 
sic split in each clone, once more because distant Fourier 
components of P^ do not affect a given solution. If 
is so large that both clones are synchronized with each 
other, a true adiabatic solution would lead to one single 
classic split for the entire spectrum, but this solution is 
unstable. Once jj, has become smaller than the separa- 
tion between the clones, the classic co-rotation solution 
for the entire spectrum gives way to two approximately 
independent solutions, one for each clone. 

Most interesting in the SN context is a triple-crossed 
spectrum. For NH we have a single unstable crossing 
flanked by two stable ones, resulting in a single swap 
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(Fig. [2]) that can be understood, using the appropriate 
in Eq. For schematic spectra one finds explicitly 
an exponential decrease of Kq with decreasing spacing of 
the crossings. For IH the crossings are unstable-stable- 
unstable and generically provide two swaps. However, 
they can merge to a single swap if the spacings between 
the crossings are too small, i.e. the triple crossing acts 
as an effective single one. The dynamical formation of 
the swaps in the triple-crossed case can be very different 
depending on the spectrum. The location and widths 
of the swaps are fixed dynamically by coupled-pendulum 
solutions together with the corresponding adiabatic in- 
variants. 

Multi-angle effects. — SN neutrinos travel on different 
trajectories and every mode is described by its direction 
of motion, in addition to w. Multi-angle effects can cause 
kinematical decoherence among angular modes [Tst and 
smear the spectral features. However, in spherically sym- 
metric cases the collective motions seem rather robust 
against multi-angle effects and a single-angle treatment 
often seems to capture the main effect 

The swap factors from multi-angle simulations of our 
examples are shown as dashed lines in Figs. [2] and [31 
The right split in Fig. [3] is robust whereas the left one 
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smears out as previously noted [20|, l2l| . In Fig. [2] we 
find significant decoherence across the spectrum. The IH 
case requires a large number of modes to reach apparent 
convergence (we used 200 frequency and 15000 angular 
modes). These issues require a dedicated investigation 
that is beyond the scope of our work. It seems clear, 
however, that some of the splits survive multi-angle de- 
coherence as sharp spectral features. 

Concluding remarks. — Spectral swaps with concomi- 
tant splits are generic for dense neutrino fluxes. They 
can appear in the v and v channels for both mass hierar- 
chies. We have provided a qualitative explanation for the 
main features of the multi-split phenomenon. Spectral 
swaps are nucleated by unstable crossings of the initial 
flavor spectra. The number and locations of observa- 
tionally relevant splits depends on the properties of the 
initial spectra and the effective mixing angle, opening the 
possibility for observing interesting time-dependent fea- 
tures. Notice that the SN flavor spectra strongly differ 
between the prompt deleptonization burst, the accretion 
phase, and the cooling phase. Multi-angle effects smear 
some of the splits, once more opening the possibility for 
interesting time-dependent features. It seems that the 
physics of SN neutrino conversion is much richer than it 
was thought before and further studies may uncover a 
number of qualitatively new phenomena. 
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